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^ ■ Abstract 
O 

The invariance for the equation of fast diffusion in the 2D coordinate space has been 
proved, and its reduction to the ID (with respect to the spatial variable) analog is demon- 
strated. On the basis of these results, new exact multi-dimensional solutions, which are 
A ■ dependent on arbitrary harmonic functions, are constructed. As a result, new exact solu- 
tions of the well-known Liouville equation - the steady-state analog for the fast diffusion 
equation with the linear source - have been obtained. Some generalizations for the sys- 
tems of quasi-linear parabolic equations, as well as systems of elliptic equations with 
Poisson interaction, which are applied in the theory of semiconductors, are considered. 

^ ■ Introduction 

The paper investigates the equation of fast diffusion 

u t = Aliau, u — u(x,y,t), (1) 

in the 2D coordinate space, which is characteristic of many applied problems, for example, 
in the description of spreading of superfine monomolecular layers of liquid under the 
influence of Van der Waals forces [1]. It arises in modeling of diffusion phenomena in 
semiconductors, polymers, etc. It is known [2] that eq. (1) is special from the viewpoint 
of the group theory since it assumes finite-dimensional algebra of point symmetries. This 
means that eq. (1) in the 2D coordinate space possesses an infinite stock of invariant 
solutions [3]. In the literature, the relation (1) is sometimes called the Ricci equation. 

This paper considers some properties of (1). On the basis of these properties some 
new - exact and multidimensional - solutions of this equation, which are dependent on 
arbitrary harmonic functions, are constructed. 
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1 Invariance of the equation of fast diffusion 



Let invariance be understood as nonvariability of the form of this equation under the 
effect of any transformation. 

Defenition[4]. Harmonic functions £(x,y) and r)(x,y) are called conjugate in the 
simple connected domain D if the function F(z) = £(x,y) + irj(x,y) is some analytical 
function of the argument z = x + iy in the domain D. 

Conjugate harmonic functions are related by Cauchy-Riemann equations 

d£ drj <9£ drj 
dx dy' dy dx' 

define one another everywhere in D with the to an additive constant and, consequently, 
possess the following properties 

(V£,Vr/)=0, |V£| 2 = IV^I 2 . (2) 

The harmonic polynomials 

V) = (x 2 + y 2 ) 2 cos(ny?), rj(x, y) = (x 2 + y 2 j 2 sin(ny?), 



where 



ip = arccos —=== , or <p — arcsm —=== , n e JM, 
\ \Jx z + y 2 J \ \Jx l + y 2 ) 



are the simplest examples of the conjugate harmonic functions. 

Lemma 1. // the function r}(x,y) is harmonic then the function In rj 2 , + rj 2 
harmonic. 



is also 



Proof. For the purpose of comfort let us introduce the denotation ip(x, y) = rj 2 , + rj 2 . 

A\ni;= [iptyxx + M-tyl + tf)]^- 2 . (3) 
By consecutive computing necessary partial derivatives, we obtain 

Ipxx + ^yy = 2 {nix + Vty) + ^Vly + 2% (Vxx + Vyy) x + (Vxx + %y) y , 
^1 + ^1 = %VxVyVxy (Vxx + Vyy) + ±vl (vl x + vly) + ^ (v^y + Vly) ■ 

Since due to the harmonic character of r) xx = —r} yy , the latter relations will, respectively, 
assume the following form: 

^xx + ^ yy = 4: (rg x + rg y ) , 
^ + ipl = 4 (rj 2 x + $ {nix + Vly) ■ 



Hence, expression (3) turns identically zero. □ 

Now let us formulate one of the main results of the work. 

Theorem 1. If £ = ^(x,y),i] = r](x,y) are conjugate harmonic functions then the 
equation of fast diffusion (1) is invariant with respect to the transformation 

u(x,y,t) = p(x,y)v(£,r],t), (4) 

where p(x,y) = |V?7| 2 , i.e. under the effect of (4) it transforms into itself 

Vt = Az v \nv. (5) 



Proof. After substitution of relation (4) into eq. (1) and simple computations, on 
account of the properties of conjugate harmonic functions (2), we obtain 

pv t = A xy In p + pA (v In v. 

Hence, due to Lemma 1, relation (5) immediately follows. □ 

Therefore, expression (4) is a form of autotransformation for the equation of fast 
diffusion (1) and, consequently, it is the formula of branching solutions. 

Example 1. Consider a exact solution of eq. (5) of the form [5, 6] 

th(t) 



By applying transformation (4) with arbitrarily chosen harmonic functions to this solution 
let us construct a few new exact nonautomodel explicit solutions of eq. (1), which are 
anisotropic with respect to spatial variables 

_ o (cth 2 (a;)+tg%))th(t) 
U{X ' V ' >~ l + cth 2 (:r)tg 2 (y)th 2 (t)' 

tg 2 (x)+th 2 (y))th(t) 



u(x, y, t) = 18 



u(x,y,t) = 18 



u (^y^)- 2 1+tg 2 {x)th 2 {y)th 2 {t y 

(tg 2 {x 3 - 3xy 2 ) + th 2 {3x 2 y - y 3 )) (x 2 + y 2 f th(t) 
1 + tg 2 (x 3 - 3x^ 2 ) th 2 {3x 2 y - y 3 ) th 2 {t) ' 

(tg 2 (exp(3x)a(y)) + th 2 (exp(3x)%))) exp(x)th(t) 



1 + tg 2 (exp(3x)a(y)) th 2 (exp(3x)6(y)) th 2 (t) 
where the following denotations are used 

a(y) = cos 3 (y) - 3cos(y) sm 2 (y) : b(y) = 3cos 2 (y) sin(y) - sin 3 (|/). 
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The solutions obtained are significant in virtue of the fact that for t — > +00 they are 
stabilized in the form of steady-state solutions. 

A result similar to that of Theorem 1 can be extended onto some class of systems of 
quasilinear parabolic equations. Hence, the following theorem is valid. 

Theorem 2. Let £ = £(x,y), rj = rj(x,y) are conjugate harmonic functions, and 
fi(ui, u 2 , ■ ■ ■ , u m ),i = 1,2, ... ,n,m < n, are homogeneous, the degree of homogeneity 
being one, i.e. 

fi{\Ui, \u 2 , • • • , \u m ) = \fi(ui, U 2 ,..., Um). (6) 

The n the system of n equations 

Ou ■ 

= A xy \nui + /i(ui,U2,...,ii m ) (7) 
is invariant with respect to the transformations 

Ui(x,y,t) = p(x,y)vi(£,ri,t), (8) 
where p = |V?7| 2 , i.e. it transforms into itself 

= Afrlnvi + fi(v u v 2 , . . . , v m ). (9) 



Proof. By substituting functions (8) into eqs. (7) and taking account of the conditions 
of homogeneity (6), we obtain the system 

dv ■ 

p-j± = A xy (In p + In Vi) + pfi(v 1 ,v 2 ,...,v m ), (10) 

furthermore, when using the properties of the conjugate harmonic functions (2), it can 
readily be shown that 

A xy \nvi = pA^lnVi. 

Therefore, from (10), due to Lemma 1, it is possible to obtain the system of equations 
(9). □ 



2 Reduction of the equation of fast diffusion to its 
one-dimensional (with respect to the spatial vari- 
able) analog 

Theorem 3. Let r](x, y) be an arbitrary harmonic function, which is different from the 
constant one. Then with the aid of the transformation 

u(x,y,t) = (pZ + rff)v(ri,t), (11) 
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the equation of fast diffusion (1) can be reduced to its one dimensional (with respect to 
the spatial variable rj) analog 

dv _ d 2 \nv 

dt drf ' 1 ' 

Proof. After substituting expression (11) into eq. (1), we obtain the equality 

dv 

ii(x,y)— = A xy (In ifj(x,y) + In v) , (13) 
furthermore, it can easily be shown that 

d 2 hif 



Axyhiv = ip(x,y)- 



drj 2 



where ip(x,y) — rj 2 + r] 2 . Therefore, from equality (13), due to Lemma 1, we obtain eq. 
(12). □ 

Consequently, by integrating eq. (12) and using formula (11), it is possible to construct 
a class of exact solutions of fast diffusion equations (1), which depend on an arbitrary 
harmonic function. 

Example 2. Let rj(x,y) be a 4th degree harmonic polynomial, i.e. n(x,y) = x 4 — 
Qx 2 y 2 + y 4 , or rj(x, y) = 4 (x 3 y — xy 3 ) . Hence r] 2 + r] 2 — 16 (x 2 + y 2 ) 3 . In this case, the 
transformation 

u(x, y,t) = 16 (x 2 + y 2 ) v(7],t), 

gives an anisotropic (with respect to spatial variables x and y) solution of eq. (1), fur- 
thermore, the function v(r],t) can be determined from the relation (12). Since with the 
aid of the substitution v = W 1 eq. (12) can be reduced to the equation with quadratic 
nonlinearities 

w t = ww m - (w v ) 2 , 

it is possible, in addition, to construct some exact solutions of eq. (12) on the basis of 
the approach described in [7]. 



1 ^k 2 + A; 2 2 sh(At) 

v(r},t) 



^ ki cos(r?) + k 2 sm(i]) + \jk 2 + k 2 ch.(\t) 

1 \Jk 2 + k\ cos(At) 

^ k\ cos(r]) + ki sin(?7) — Jk\ + k 2 sin(At) 



v(rj,t) 



1 \/k 2 - k\ cos(At) 



^ hchir]) + k 2 sh(r]) + y/k% - k\ sin(At) ' 
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1 sjkl - Afah(At) 



A fcidi^) + A; 2 sh(?7) - - fcfch(At) 
Here A ^ 0, fci > ^ are arbitrary numerical parameters. 

Remark. The proof of an analog of Theorem 3 for the system of quasilinear parabolic 
equations (7) has been described in the author's paper [8]. 

It is known rather well [9, .84] that the Laplace equation in the 2D coordinate space 
is invariant with respect to the conformal transformation of independent variables 

dx dy dx dy 



dq 2 dqi ' dqi dq 2 ' 



(14) 



dx dy 

furthermore, — — , — — do not simultaneously turn zero. The system of parabolic coordi- 
oqi oq 2 

nates 

x = ql - qf, y = 2q ± q 2 . 

is an example of such transformation. Now, let us apply the transformation (14) to eq. 
(1) under scrutiny. The latter may be rewritten as follows: 

u t = /(gi,g2)A 9l q 2 lnM, u = u(q 1 ,q 2 ,t). (15) 

Here the denotation f(qi,q 2 ) = x 2 qi + x 2 q2 is used. Since due to (14) the function x = 
x(qi,q 2 ) is harmonic, Lemma 1 implies that 

A qiq2 \nf( qi ,q 2 )=0. (16) 

Consequently, the following theorem is valid. 

Theorem 4. If relation (16) holds then with the aid of the transformation 

u(qi,q 2 ,t) = — - \f 2 qi + f 2 } v(r},t), rj = In f( qi , q 2 ), (17) 

the nonhomogeneous equation of fast diffusion (15) can be reduced to the homogenous and 
one- dimensional (with respect to the spatial variable n) equation (12). 

The proof of this theorem is similar to that of Theorem 3. To the end of proving it is 
necessary to put r)(qi,q 2 ) = \nf(qi,q 2 ) in (11) written in terms of variables q\ and q 2 . 

Example 3. The equation of nonlinear diffusion 

Ut = exp(x 2 — y 2 ) A In u 

or w t = exp(x 2 — y 2 — w)Aw, 

where w(x,y,t) = \nv(x,y,t), has the following exact solution, which is asymmetric with 
respect to the spatial variables 

u(x, y, t) = A(x 2 + y 2 ) exp(x 2 — y 2 )v(rj, t), rj = x 2 — y 2 . 

Furthermore, the function v(r],t) satisfies eq. (12), some exact solutions of which have 
been given in Example 2. 
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3 Exact solutions of Liouville equation in the 2D co- 
ordinate space 



It can easily be shown that principal results of the sections 1 and 2 can be directly 
generalized onto the equation of fast diffusion with a linear source (sink) 



u t = A In it — Ait, 



(18) 



where A G 9? \ {0}. Since formulas (4), (11), (17) are obviously independent of time, 
Theorems 1, 3 and 4 remain valid for the steady-state (elliptical) eq. (18) 

A In u — Ait, 

which by the substitution u = exp(\w) may be reduced to the well-known Liouville 
equation: 

Aw = exp(Xw). (19) 

So, Theorem 1 and Lemma 1 together imply the following 

Proposition 1. The Liouville equation (19) is invariant with respect to the transfor- 
mation 

1 

w(x,y) = v(£,ri) + -\n\p(x,y)\, (20) 

where p — |V^| 2 , and £(x,y),r)(x,y) are conjugate harmonic functions. 

By applying transformation (11) to eq. (19) it is possible to obtain an ordinary 
differential equation (ODE), which can be easily integrated, and make sure that the 
following proposition is valid. 

Proposition 2. The Liouville equation (19) in the 2D coordinate space has exact 
solutions of the form 



w(x,y) = - In 



2A 2 



w(x,y) 



A 



In 



A 
2A 2 



A 



(vl + vl) sec 2 {Ar]) 
(pi + sech 2 (A V ) 



where rj(x,y) is an arbitrary harmonic function, which is different from the constant one, 
A,X E K\{0}. 

Similarly, from Theorem 4 it follows that for a steady-state nonhomogeneous equation 
(15) having a linear adjunct the following proposition is valid. 



Proposition 3. The nonhomogeneous Liouville equation 

Aw = r](x, y) exp(Aio), 



(21) 
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where r)(x, y) is an arbitrary harmonic function different from the constant one, has the 
exact solution 

w(x, y) = — In 



A 



3vl + vl 



A ?7 3 



Note that exact solutions of the nonhomogeneous Liouville equation (21), when r](x, y) 
is a holomorphic function of special form, can be found in [10]. 

As a case of generalization of the results obtained in this section onto some different 
systems of equations consider the steady-state mathematical model of charge transfer, 
which implies Poisson interactions, known from the theory of semiconductors. 

Example 4. The elliptical system of equations of the form 

Au = exp(w) + A\V\ 2 , Av = exp(v) - B\V\ 2 , 

A = exp(i;) - exp(it), (22) 

has a exact solution of the form 

u(x, y) = f(rj) + In \-q 2 x + r/j| , v(x, y) = ^(r/) + In \rf x + 7/j| , (x, y) = <p(rj), 

where rj(x, y) is an arbitrary harmonic function, which is different from the constant one; 
A, B G 3? \ {0}. Furthermore, the functions /, if), <p can be defined from the system of ODE 

/" = exp(/) + V, r = exp(^) - B/, ip" = exp(^) - exp(/). 



In the case when A = —B, system (22) assumes a particular exact solution u(x, y) 

v{x,y) = f(rj) + In 
linear 2nd order ODE 



, (x,y) = r](x,y), furthermore, the function f{rj) satisfies the 
f" = f + A. 
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